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Abstract. 

Radiative transfer equilibrium models of nearby interstellar matter (ISM) are required 
to determine the boundary conditions of the heliosphere from astronomical observations 
of nearby stars. These models are also constrained by data on the ISM inside of the so- 
lar system, including pickup ion, anomalous cosmic rays, and in situ He° data. The two 
best ISM models give semi-empirical filtration factors for H (0.41 - 0.52), He (0.86-1.16), 
N (0.74-1.11), O (0.55-0.77), Ne (1.41-2.80), and Ar (0.61-1.76) when observational un- 
certainties are included. Uncertainties in the Ne filtration factor may result from poorly 
known interstellar abundances. These models predict the characteristics of the ISM out- 
side of the solar system: n(H°)=0.20-0.21 cm" 3 , n(e~)=0.10 cm" 3 , H+/H=0. 29-0.30, 
and He + /He=0.47-0.51. However, if the isotropic 2 kHz emission observed by Voyager 
(Kurth & Gurnett 2003) is formed in the surrounding ISM, an alternate model (Model 
25) is indicated. The weakly polarized starlight of nearby stars suggests that the local 
galactic magnetic field is parallel to the galactic plane, and the strongest polarization 
is towards the upstream direction of the ISM flow, and also (coincidently) near the eclip- 
tic plane. Observations of nearby ISM, the radiative transfer models, and historical 10 Be 
records provide information on past variations in the galactic environment of the Sun. 



1. Introduction 

The boundary conditions of the heliosphere are domi- 
nated by the ionization, density, temperature, and mag- 
netic field of the surrounding interstellar cloud. The first 
use of anomalous cosmic ray (ACR) data to study in- 
terstellar ionization compared ACRs and the Mg°/Mg + 
ratio towards Sirius, relying on sightline-averaged ion- 
ization rates (Frisch 1994). Full radiative transfer equi- 
librium models of interstellar gas within ~3 pc are now 
available, predicting ISM properties at both the helio- 
sphere boundary and averaged over sightlines towards 
nearby stars (Slavin & Frisch 2002, SF02, Frisch & Slavin 
2003, FS03). Such models provide a tool for evaluating 
the boundary conditions of the heliosphere when con- 
strained with interstellar matter (ISM) data inside and 
outside of the heliosphere. However, the loss of interstel- 
lar neutrals to charge exchange with interstellar ions in 
the heliosheath regions ( "filtration" , Ripken & Fahr 1983, 
Izmodenov et al. 1999, Mueller & Zank 2002, Cummings 
& Stone 2002, CS02) must be considered when interpret- 
ing the in situ data. In this paper, the radiative model 
predictions for the physical properties of the ISM at the 
solar location are used to determine semi-empirical filtra- 
tion factors, extending the discussion presented in FS03. 
Models consistent with a possible interstellar origin for 
the 2 kHz emission observed by Voyager in the outer 
heliosphere (Kurth & Gurnett 2003) are also discussed, 
along with evidence for a nearby interstellar magnetic 
field and time-variability of the boundary conditions of 
the heliosphere. 



2. Why Radiative Transfer Models? 

The ionization of the ISM at the solar location and 
boundary conditions of the heliosphere are highly sensi- 



Copyright 2008 by the American Geophysical Union. 
0148-0227/08/$9.00 



1 



X - 2 



FRISCH : HELIOSPHERE-ISM INTERACTION 



tive to the interstellar cloud properties integrated from 
the solar location (SoL) to the cloud surface. Since the H° 
column density to the local cloud surface is log iV(H°)<18 
cm~ 2 , and photons which ionize H (>13.6 eV) and He 
(>24.6 eV) are attenuated to 1/e by log iV(H°)~17.2 
cm~ 2 and log iV(He )~17.7 cm -2 respectively, the H° 
and He° densities at the SoL are sensitive to the radia- 
tion attenuation inside the cloud and to cloud geometry. 
The measurements which best distinguish the properties 
of the surrounding cloud are those which sample par- 
tially ionized elements with first ionization potentials in 
the range 13-25 eV (e.g. H, He, O, N, Ne, and Ar), which 
also constitute the parent population of pickup ions and 
anomalous cosmic rays. 

Radiative transfer (RT) models of ISM near the Sun 
(<3 pc) are required to predict the ionization of the 
ISM at the SoL. RT models have been constructed us- 
ing measurements of the interstellar radiation field in the 
interval 300-3000 A, a modeled conductive interface be- 
tween nearby ISM and the hot plasma (~10 6 K) inte- 
rior of the Local Hot Bubble, and the cloud composi- 
tion determined from observations of ISM both inside 
and outside of the solar system (see SF02 and FS03). 
Free parameters in these models include the H° column 
density towards the cloud surface (iV(H°), which is not 
directly measured for the sightline used to constrain the 
models, e CMa), the soft X-ray emission from the Local 
Bubble hot plasma (characterized by plasma temperature 
Th), the interstellar magnetic field strength (B, which 
regulates conduction at the cloud boundary and the ex- 
treme ultraviolet radiation field, Slavin 1989), and cloud 
density (nH=n(H°)+n(H + )). Sightline-averaged values 
of the relative ionizations of H and He in the ISM are 
found from observations of the He° ionization edge (504 
A) towards nearby white dwarf stars (40-100 pc), which 
are strong sources of radiation in this extreme ultravio- 
let (EUV) wavelength interval (e.g. Kimble et al. 1993, 
Vallerga 1996). These data generally sample sightlines 
with several blended clouds, and show a wide variation 
in ionization (JV(H°)/iV(He o )~10-15). The equilibrium 
models fitting the high He ionization seen in the EUV 
data require an excess of radiation in the 500 A region, 
when compared to the soft X-ray data, and the variation 
of H°/He° with cloud depth requires radiative transfer 
models (e.g. Cheng & Bruhweiler 1990, Kimble et al. 
1993, Vallerga 1996, SF02). 

The initial set of 25 models was pared to seven "best" 
models (Frisch and Slavin 2003, FS03) by introduc- 
ing recent astronomical data indicating gas-phase in- 
terstellar abundances of O/H~400 ppm in long (>700 
pc) interstellar sightlines through predominantly diffuse 
clouds (Andre et al. 2002), since O and H ioniza- 
tions are closely coupled by charge exchange. Data on 
Mg°/Mg+ and C + */C+ towards e CMa (Cry & Jenkins 
2001) then indicated that Models 2, 8, and 18 provide 
the best match to astronomical data, and Model 11 is 
within uncertainties (FS03). FS03 compared these mod- 
els with He data inside of the solar system, and con- 
cluded Models 2 and 8 are the best models if He is not 
filtered in the heliosheath regions. (Ionization and den- 
sity predictions for the full 25 models are available at 
http:/ /astro. uchicago.edu/~frisch/SlavinFrisch2002_tables.) 
The principal failures of these two RT models are that 
predicted Ne° densities are ~50% of Ne° densities in- 
ferred from PUIs and ACRs, and cloud temperatures are 
~30% larger than values inferred from the He° data. The 
filtration predicted by these models is examined below. 
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3. Filtration Factors 

Gloeckler & Geiss (2001, GG01) made an empirical es- 
timate of hydrogen filtration, Fh = [i%Ts(H)/nTs(He)]/(j), 
by comparing H and He densities at the termination 
shock with the interstellar values, 0=nis,Sun(H°)/nis,Sun(He°), 
assuming negligible He filtration. They assumed 
0=11.25 based on white dwarf data which show 
JV(H°)/iV(He )~11.25, yielding Fn= 0.54. GG01 also 
derived an O filtration factor Fo= 0.62, using H° and 0° 
column densities towards the stars Procyon and Capella. 
(Volume densities are n, in units cm - , and column den- 
sities for element X are N(X), in units cm -2 ). The SF02 
models show iV(H°)/iV(He )=8.8-13.6, consistent with 
the large variations found for N(H°)/N(He°) towards 
nearby stars, and predict values 0=12. 1-14.7 at the SoL. 
The relative H and He ionizations at the SoL depend 
on the model parameters (4>=4>(nH,N(H°),B,Tn)), since 
H-ionizing photons are attenuated more strongly than 
He-ionizing photons. 

Filtration factors can be determined semi-empirically 
from comparisons between ionization predictions for the 
SoL (Table 7 in SF02) and in situ data on the byprod- 
ucts of ISM interactions with the solar system. Neu- 
tral densities at the termination shock are determined 
from pickup ion (GG03), anomalous cosmic ray (CS02), 
and direct He measurements (Witte et al. 2003). 
These data give n T s(He°)=0.0145±0.0015 cm -3 and 
n T s(H°)=0.095±0.01 cm -3 . The H and He filtration fac- 
tors determined from the 25 models are shown in Fig. la, 
along with theoretical predictions of H and He filtration 
factors. The seven models with O/H~400 ppm are plot- 
ted as filled circles, or triangles for models 2, 8, 18 which 
provide the best match to astronomical data. Among 
the three models best matching astronomical data, 2 and 
8 are most consistent with predictions of negligible He 
filtration. The Fn e predictions for the upstream direc- 
tion are F He =0.98-0.99 (Mueller and Zank, 2002) and 
F Hc =0.94 (CS02). The hydrogen filtration factors are 
more uncertain, F H =0.4-0.5 (Mueller & Zank 2002, Rip- 
ken & Fahr 1983, Izmodenov et al. 1999, Gloeckler & 
Geiss 2003). 

Empirical filtration factors for Ar, O, and N are also 
shown in Fig. 1, and in all cases Models 2 and 8 re- 
sults are consistent with the predicted range of values 
(when uncertainties are included). Model 18 is not con- 
sistent with theoretical calculations suggesting Fh c ~ 1, 
but is consistent with Ar, O, and N predicted nitrations. 
Models 2 and 8 provide the best fit to the combined in- 
terstellar (Mg and C, FS03) and heliospheric data, and 
predict similar boundary conditions for the heliosphere 
(n(H°)~0.21 cm -3 , n(e")~0.10 cm -3 , see Table 1). The 
radiative transfer models are not well constrained us- 
ing heliospheric data alone. The semi-empirical filtra- 
tion factors predicted by Models 2, 8 are, respectively: 
H (0.46±0.05, 0.47±0.05), He (0.95±0.09, 1.06±0.10), 
N (0.91±0.17, 0.93±0.18), O (0.65±0.10, 0.67±0.10), Ne 
(1.75±0.34, 2.34±0.46), and Ar (1.10±0.49, 1.22±0.554, 
for Ar/H=3.16e-6, see below). The filtration factors 
based on Model 11 are within the range quoted for Model 
2, although the uncertainties differ. 

The result that Jnc >1 indicates that either the inter- 
stellar Ne abundance is poorly understood, or that a he- 
liosheath process enhances Ne° (such as charge exchange 
with He°, with a similar ionization potential). The poorly 
understood extreme ultraviolet radiation field (570-790 
A, EUV) ionizes both Ne° and He°, and since He° data 
are consistent with model predictions, uncertainties in 
the EUV radiation field do not appear as the likely expla- 
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nation. The interstellar Ar abundance is uncertain (Sofia 
& Jenkins, 1998), and the abundance Ar/H=3.16e-6 is 
adopted here. 

As a consistency check, the ratio of the PUI n(0)/n(N) 
to the ISM N(0)/N(N) values have been calculated, and 
compared to the ratio Fo/Fn. For the ISM, the ra- 
tio AT(0°)/JV(N°) towards e CMa is used (SF02, Gry 
& Jenkins 2001) (although this is not equal in de- 
tail to n(0°)/n(N°) at the solar location because N° 
has a higher photoionization cross section than 0°). 
The results are: (ra(0)/n(N) pui )/(iV(0)/iV(N) iBm ) = 
0.70±0.28, compared to F o /F N =0.72±0.20, as expected. 

4. 2 kHz Radio Emission 

Gurnett et al. (1998, 2003) have detected a relatively 
isotropic radio-emission signal at 2 kHz in the outer he- 
liosphere with the plasma wave instruments on board 
the two Voyager spacecraft; a possible explanation for 
the isotropy is emission at the plasma frequency of the 
surrounding interstellar gas providing this emission can 
enter the heliosphere (e.g. Izmodenov et al. 1999). If 
this emission has an interstellar origin, n(e~)~0.05 cm~ 
and the closest match to this electron density among 
the SF02 models would be Model 25 (n(e~)~0.0607 
cm~ 3 ) and Model 20 (n(e~)~0.0618 cm -3 ). Model 
25 predicts: n(H°)=0.21 cm -3 , n(He°)=0.017 cm" 3 , 
H+/H=0.194, He+/He=0.341 and T=5,120 K. Model 
20 predicts: n(H°)=0.23 cm -3 , n(He°)=0.018 cm' 3 , 
H+/H=0.183, He+/He=0.346 and T=7,200 K. The rel- 
atively low ionization levels found by these models re- 
sults from the absence of a conductive interface on the 
cloud, which is characterized by the magnetic field pa- 
rameter B=0. Model 25 yields predicted filtration fac- 
tors of: F H =0.45±0.05, F Hc =0.85±0.08, F N =0.55±0.11, 
F o =0.39±0.06, F Ne =1.16±0.23, and F Ar =0.83±0.37. 
GG03 have found that Model 25 offers a self-consistent 
solution when heliosphere models of filtration are com- 
pared to the PUI data. However, if Model 25 is correct 
further study is required to understand the local ISM to- 
wards e CMa, since the Model 25 predictions for interstel- 
lar abundances do not match the local ISM observations 
towards e CMa. 

5. Interstellar Magnetic Field 

The interstellar magnetic field (Bis) immediately out- 
side of the heliosphere has not been directly measured. 
However, observations of weakly polarized light from 
nearby (<50 pc) stars show a patch of interstellar dust 
centered on galactic coordinates Z=350°-^ Z=20°, b=- 
40°^ b=-5°, with extinction A Y ~0.01 mag (or iV(H)~2 
x 10 19 cm -2 ) and galactic field directed towards / ~70° 
(Tinbergen 1982, Frisch 1990). Fig. 2 shows the location 
of the stars (in galactic coordinates), and the polarization 
vectors, in the region where the weak polarization is most 
evident, with the plane of the ecliptic superimposed. Also 
plotted are the radio emission sources detected by the 
Voyager spacecraft which are consistent with an origin 
near the heliosphere nose and follow roughly the galactic 
plane (Kurth & Gurnett 2003). The weak polarization is 
most evident near the plane of the ecliptic and coincides 
with the upstream direction of the cluster of interstel- 
lar clouds flowing past the Sun (which is shown in Fig. 
2 plotted in the Local Standard of Rest, after correc- 
tion for solar motion). These polarizing dust grains are 
predominantly within 5 pc of the Sun, or closer (Frisch 
1990). Although interstellar grain alignment mechanisms 
are poorly understood, popular grain alignment models 
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yield polarization vectors parallel to the interstellar mag- 
netic field (Lazarian & Cho 2002). The same classical in- 
terstellar dust grains which polarize optical radiation also 
pile up in the heliosheath regions as grains are deflected 
around the heliosphere (grain radius ~0.1 /jm, Frisch et 
al. 1999). Rand & Lyne (1994) compared pulsar disper- 
sion and rotation measures towards pulsars spaced over 
several kpc, to derive the strength of the ordered compo- 
nent of the magnetic field near the Sun, B=1.4±0.2 jj,G, 
directed towards galactic longitude I = 88° ±5°; this field 
is consistent with the field traced by Tinbergen's data. 
Once the production mechanisms for the grain alignment 
and for the radio emission sources are understood, it may 
produce a deeper understanding of the asymmetries in 
the heliosphere nose region due to the ~60° tilt between 
the ecliptic plane and interstellar magnetic field. 

6. Variability in the Boundary Conditions 
of the Heliosphere. 

Simple geometrical assumptions indicate the Sun left 
the Local Bubble interior and entered the surrounding 
interstellar cloud within the past several thousand years 
(Frisch 1994). More precise estimates of the variability 
of the boundary conditions of the heliosphere result from 
combining densities from the best RT models with obser- 
vations of nearby stars. 

Models 2 and 8 give a local ISM density of n(H°)~0.21 
cm" 3 , which when combined with the column density to 
the cloud surface (JV(H°)=6.5 10 17 cm" 2 ) give a distance 
to the cloud surface of ~0.97 pc. This distance is tra- 
versed in ~37,000 years by a cloud at the LIC velocity. 
Since the local ISM towards e CMa and a CMa (Sirius) 
consists of two cloudlets (Gry & Jenkins 2001, Lallement 
et al. 1994), and the second blue-shifted cloud has a 
higher velocity relative to the Sun than the LIC, the time 
at which the Sun encountered the LIC must be shorter 
than 37,000 years. Frisch (1997) attributed spikes in the 
10 Be ice core record 33 kyr and 60 kyr ago to encoun- 
ters with ~0.1 pc wide magnetic structures, comoving 
in the local ISM, which confine low energy galactic cos- 
mic rays. These structures would now be ~0.87 pc, and 
~1.6 pc distant from the Sun in the downstream direc- 
tion (e.g. towards e CMa, a CMa). The most recent 10 Be 
spike may be consistent with the solar entry into the local 
ISM; alternatively it may correlate with the Sun's entry 
into the LIC while the earlier event may signal entry into 
the blue-shifted cloud seen towards CMa. Prior to that, 
the Sun was immersed in the hot (10 4 K) plasma interior 
to the Local Bubble where an extended heliosphere (he- 
liopause radius 300 au in the nose direction) is predicted 
(Frisch et al. 2003). (Civilization developed as the Sun 
emerged from the Local Bubble interior, where it resided 
for several millions of years.) 

Looking to the future, at least three interstellar clouds 
are found within 5 pc of the Sun in the upstream direc- 
tion, the Local Interstellar Cloud (LIC) where the Sun is 
now located, the "G-cloud", and the Apex Cloud (Frisch 
2003). The velocity vector of each cloud is consistent with 
the cloud being located in front of the nearest star a Cen, 
although the G-cloud is the best fit to observations (e.g. 
Lallement et al. 1995, Linsky & Wood 1996) and the LIC 
would have to have a small velocity gradient (~0.7 km/s) 
towards the star. If the G-cloud is the next cloud to be 
encountered by the Sun, the higher heliocentric velocity 
(-29.1 km/s, factor of ~1.1 ) and neutral density (>5 
cm" 3 , factor of >24) than for the LIC indicates a dra- 
matic heliosphere shrinkage if the other properties (e.g. 
ionization) are similar to the LIC (Frisch 2003, Zank & 
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Frisch 1999). If, in contrast, the Apex Cloud is the next 
cloud encountered, the larger heliocentric velocity (-35.1 
km s" 1 ) also indicates a smaller heliosphere, although the 
other properties of this cloud are unknown. Variations of 
the boundary conditions of the heliosphere are thus ex- 
pected on timescales of less than 37,000-45,000 years if 
either the G-cloud or Apex Cloud is foreground to a Cen. 

7. Discussion 

The primary results of this paper are estimates of fil- 
tration factors for H, He, N, O, Ne, and Ar, based on 
radiative transfer models of the local ISM from SF02, 
and FS03. Two radiative transfer equilibrium models of 
ISM within ~3 pc of the Sun (models 2 and 8) yield 
predictions which match observations of ISM inside (He, 
PUIs, ACRs) and outside (gas within ~3 pc towards e 
CMa) of the solar system, for an assumed interstellar 
gas-phase O abundance of O/H=400 ppm (FS03). These 
models yield semi-empirical predictions for He, H, O, N, 
Ar filtration factors which are consistent with theoretical 
predictions, within uncertainties, and predict interstellar 
densities n(H°)~0.20 cm" 3 , n(e")=0.10 cm" 3 at the so- 
lar location. However, there are unexplained differences 
between the predictions of Models 2 and 8 and observed 
pickup ion Ne densities and the cloud temperature. 

If the astronomical constraints on the radiative trans- 
fer models are relaxed, alternate models become plausi- 
ble. For example, if the 2 kHz radio emission observed 
by Voyager 1, 2 is from the surrounding ISM, the low 
electron density of Model 25 (n(e")~0.06 cm" 3 ) is a 
better match although predicted interstellar abundances 
are unrealistic. For Model 25 to be correct the radia- 
tion field intensity versus cloud geometry requires fur- 
ther study. The simple interstellar cloud geometry as- 
sumed for the underlying radiative transfer models may 
also prove inadequate. If the astronomical constraints 
prove incorrect (e.g. the assumed interstellar gas-phase 
oxygen abundance), then the selection of the best mod- 
els to describe the boundary conditions of the heliosphere 
may also change. 

As better data on the ISM inside and outside of 
the heliosphere become available, yielding more accu- 
rate boundary conditions for the heliosphere, heliosphere 
modeling will predict more reliable filtration factors for 
comparison with semi-empirical filtration factors utilizing 
interstellar modeling. 

The strength of the interstellar magnetic field outside 
of the heliosphere remains uncertain, but both low fre- 
quency radio emission events and the polarization of the 
light of nearby stars suggest the field orientation is par- 
allel to the galactic plane and not the ecliptic plane, in- 
troducing a north/south heliosphere asymmetry. 

The second important conclusion of this paper is that 
the boundary conditions of the heliosphere change with 
time, and on timescales of <10 4 years, indicating past 
and future variations in heliosphere dimensions and prop- 
erties. 
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Figure 1. Empirical filtration factors for H (a), Ar 
(b), N (c), and O (d) are plotted against FH e , as de- 
rived from RT models (SF02, Table 7) compared to PUI, 
ACR, and He data (including data uncertainties). Filled 
circles and triangles are ISM models with O/H~400, and 
the labeled triangles are nos. 2, 8, 18 which match local 
ISM Mg+/Mg° and C+/C+* data. The box shows Model 
25 (consistent with a low electron density) . Dashed lines 
are theoretical filtration factors. 
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Figure 2. Plot of an indicator of the nearby interstellar 
magnetic field, in galactic coordinates. The bars show 
the direction of the electric vector polarization, which is 
parallel to the interstellar magnetic field direction, for 
several nearby stars (Tinbergen 1982). The arrow shows 
the likely direction of the interstellar magnetic field near 
the Sun based on the polarization and pulsar (Rand & 
Lyne 1994) data. The curved line is the ecliptic plane. 
The observed region of maximum polarization is near the 
ecliptic plane, but is also in the upstream direction (as 
referred to the "LSR" velocity rest frame of stars close 
to the Sun, Frisch et al. 2002) of the local flow of ISM 
past the Sun (star). The heliosphere nose direction, in 
the heliocentric rest frame, is located at the box. For 
comparison, the sources determined for the ~3 kHz radio 
emissions detected by Voyager are plotted as filled boxes 
(Kurth & Gurnett 2003). 
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Table 1. ISM Physical Properties at the Hcliosphere" 

Quantity Models 

2 8 



Assumed Model Parameters 



n H (cm- 3 ) b 


0.273 


0.273 


log T h (K) 


6.0 


6.1 


B (jiG) 


5.0 


5.0 


N H ° (10 17 cm- 2 ) c 


6.5 


6.5 


Predicted Quantities at Solar Location 


n(H°) (cm- 3 ) 


0.208 


0.202 


n(He°) (cm" 3 ) 


0.015 


0.014 


n(e) (cm -3 ) 
x(H) d 


0.098 


0.101 


0.287 


0.300 


X (He) d 


0.471 


0.511 


T (K) 


8,230 


8,480 


ri(N°)/n(He°) 


5.64e-4 


6.12e-4 


n(0°)/n(He°) 


5.35e-3 


5.76e-3 


n(Ne°)/n(He°) 


3.86e-4 


2.37e-4 


n(Ar°)/n(He°) e 


1.21e-5 


1.21e-5 


Predicted Quantities towards 


e CMa 


log iV H (cm- 2 ) 6 ' c 


18.03 


18.02 


A^(H )/Ar(He°) c 


11.63 


12.74 



11 From (Slavin & Frisch 2002) and (Frisch & Slavin 2002). 
6 n H = n(H°)+n(H+) cm-' 3 ; N H = N(B.°)+N(H+) cm- 2 . 

c ^(H ) is the H° column density between the Sun and cloud surface, etc., where the cloud is the sum of the two nearby cloudlets 
(d<3 pc) observed towards e CMa (Gry & Jenkins 2001). 

d x(H)> x(H e ) are the ionized fractions of H, He, respectively. 

e The B-star Ar abundance of 3.16 ppm is used here (from Sofia & Jenkins 1998). 



